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Abstract 

A tantalizing version of Maxwell's demon is presented which ap- 
pears to operate reversibly. A container of hard core disks is sep- 
arated into two chambers of equal volume by a membrane that 
selects which disk can penetrate depending on the disk's angle 
of incidence. It is shown that the second law of thermodynam- 
ics requires the incompressibility of microscopic dynamics or an 
appropriate energy cost for compressible microscopic dynamics. 



1 Introduction 
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The second law of thermodynamics is usually attributed to the fact that 
states of maximum disorder in a statistical system have the largest probability 
of occurrence among all possible states. The microscopic dynamics of a 
statistical system are assumed to be of minor importance. However the 
microscopic dynamics can not be arbitrary, and they must satisfy certain 
conditions for the second law to hold. In particular they must conserve phase 
space volume. Hamiltonian dynamics conserve phase space volume according 
to Liouville's theorem. However the condition of incompressibility is more 
general than the condition of Hamiltonian dynamics, and it is important 
to discuss incompressibility and the second law directly without the aid of 
Hamiltonians. 

This paper examines the relation between compressible microscopic dy- 
namics and the second law of thermodynamics using a time-reversible system 
of hard core disks and a membrane. The possibility that a Maxwell's demon 
could imitate the membrane of our system in a reversible manner is also 
explored, and it is shown that a Maxwell's demon can only approximate 
the membrane using irreversible operations. The membrane selects which 
disk can penetrate depending on the disk's angle of incidence, and creates 
a density difference between two chambers initialized to have equal density. 
The analysis of our system confirms that the second law of thermodynam- 
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ics requires the incompressibility of microscopic dynamics or an appropriate 
energy cost for compressible microscopic dynamics [1]. 

The next section describes our system and explains how the membrane 
creates a density difference between the two chambers. Section 3 estimates 
the density difference theoretically, and section 4 confirms the theoretical 
estimate using a computer simulation. Section 5 explains how the membrane 
challenges the second law of thermodynamics, and discusses the evolution 
of the system of disks in phase space. Section 6 shows that the membrane 
compresses the phase space of the disks, and section 7 demonstrates that 
if the membrane interaction obeyed incompressible dynamics, then it would 
not create a density difference (for piecewise differentiable maps). 

Section 8 reviews another system that obeys compressible dynamics: a 
pump of disks consisting of a microscopic trapdoor and a cooling mechanism. 
Section 9 analyzes the membrane as a Maxwell's demon that interacts with 
the disks based on information about the disks. The new demon uses a 
tennis racket to bounce the disks and to select which disk can penetrate the 
membrane depending on the disk's angle of incidence. 

2 Description of the Model 

Our system is shown in figure f. It consists of a box containing hard core 
disks, and a membrane that separates the box into two chambers of equal 
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volume. The hard core disks move and coUide with each other elastically. 
The membrane interacts with each incident disk according to the following 
equation, where Vi,V2 are the velocity components of the disk before the 
interaction and V/, V2 are the velocity components of the disk after the in- 
teraction. 



The labels +45, —45, and +90 are motivated by the discussion of section 9 
that views the membrane as a Maxwell's demon playing tennis with the disks. 



figure 2) reverses the x-component of its velocity when it hits the membrane 
and is not allowed to penetrate. It also says that an incident disk in the 
remaining octants penetrates the membrane by being defiected toward the 
X-axis when coming from the left, and away from the x-axis when coming 
from the right. The speed and the kinetic energy of the disk are conserved. 

Equation 1 is illustrated in Figure 2. The vertical solid line in figure 2 
denotes the membrane, and the other solid lines denote a division of the plane 
into octants. The two dashed lines inside the first and sixth octants (counting 




(1) 




otherwise 



The above equation says that an incident disk in octants 2,4,5, and 7 (refer to 
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counterclockwise) denote trajectories that are deflected and penetrate the 
membrane according to the case +45 of equation 1. The transformation of 
the velocity takes place instantaneously when the center of an incident disk 
reaches the membrane. 

Another way of examining equation 1 is to rewrite the equation as a map 
of the velocity angle (impact angle). This is shown in flgure 3 where 
is the impact velocity angle and 0' is the transformed velocity angle after 
the membrane interaction has occurred. Both and 0' range from — tt 
to TT. The map consists of eight line segments which can be interpreted 
as follows: A completely transparent membrane is a straight line of slope 
"1" passing through the origin, and a completely impenetrable membrane 
(mirror-reflecting) is a line of slope "—1" moved upwards tt radians away 
from the origin and wrapped around periodically. The eight line segments of 
flgure 3 can be interpreted as breaking a line of slope "—1" (mirror-reflecting) 
and shifting some of the broken pieces up and down. The broken pieces that 
are shown dashed in flgure 3 correspond to impact angles that penetrate the 
membrane. The broken pieces that are shown solid in flgure 3 correspond to 
impact angles that are mirror-reflected back. 

The membrane of equation 1 leads to a large density difference between 
the two chambers. The success of the membrane depends on the thermal 
motion of the disks and the impact rate of disks hitting the membrane. 
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Assuming that the velocities of the disks are distributed isotropically inside 
the container, it follows from geometrical considerations that the impact 
rate of disks hitting the membrane is a cosine of the impact angle in absolute 
value, 



The membrane exploits this cosine distribution of impact angles by allowing 
disks with "high rate" impact angle to penetrate from the right, and allowing 
disks with "low rate" impact angle to penetrate from the left. To achieve 
reversibility, the remaining impact angles ("low rate" from the right and 
"high rate" from the left) are blocked and do not penetrate the membrane. 
They are simply reflected back. Furthermore, each "high rate" angle from 
the right is rotated onto a "low rate" angle when the disk penetrates the 
membrane, and vice versa. 

The membrane interaction of equation 1 makes the membrane more per- 
meable from the right side than from the left. This leads to an excess flux 
of disks from the right side, and creates a density difference between the two 
chambers. When the density difference reaches an equilibrium value (which 
is calculated in the next section), the fluxes of disks between the two sides 
of the membrane become equal. 



impact-rate oc 



cos 
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3 Estimate of the Density Difference 

At equilibrium the fluxes of disks from the left and from the right side of the 
membrane are equal to each other. In other words if Nl is the normalized 
density in the left chamber, and P(^l^ji^ is the probability of an individual 
disk to penetrate the membrane coming from the left, we require that 

Nl P(l^r) = Nr P^n^L) . (3) 

We can estimate P{^l^b) using the fact that among all the disks that strike 
the membrane from the left only those with trajectories in the third and 
sixth octants of flgure 2 are allowed to penetrate. In particular, these disks 
have impact angles in the intervals (7r/4,7r/2) and (7r,57r/4). In addition, 
we assume that the probability of an individual disk to strike the membrane 
varies as the cosine of the impact angle and is independent of the density in 
each chamber. Although the total impact rate depends linearly on the density 
of the disks, the probability of an individual disk to reach the membrane is 
independent of the density to a flrst approximation. Thus we get, 

P(L^K) ~ 2 C cos Ode ~ 2 C 0.3 , (4) 

for some normalization constant C ; and similarly, 

P(R^L) ~ 2 C r'\osede ~ 2 C 0.7 , (5) 

which gives 

Nl ~ 0.7 and Nr ~ 0.3 . (6) 
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In other words the system of membrane and disks reaches equihbrium when 
the fluxes of disks from the left and from the right side of the membrane are 
balanced, and this happens when the normalized density is approximately 
0.7 in the left chamber and approximately 0.3 in the right chamber. 

4 Simulation Results 

To check the theoretical results of the previous section, a two-dimensional 
system of hard core disks with a membrane has been simulated. The com- 
puter program used in these simulations is the same program as the one 
described in detail in reference [2] with a few modiflcations to simulate the 
membrane. In particular, when the center of an incident disk reaches the 
membrane, the disk's velocity is transformed according to equation 1. 

In the experiments reported below forty disks are used. The size of each 
chamber is 24.3 X 10"^"^ cm^ (equal size chambers), and the disk radius is 
3 X 10~®cm. These numbers give a mean free path of the order of 10~^cm 
which is approximately the length of each chamber. The average speed of 
each disk is 3.56 X lO'^cm/ sec. 

In flgure 4 we can see that the time average of the number of disks in the 
left chamber increases from an initial value of 0.5 (normalized) to a value of 
0.7 as a result of the membrane interaction. The smooth curve of flgure 4 
plots the cumulative time average of the number of disks in the left chamber. 
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which approaches a steady value as the averaging time increases. The noisy 
curve of figure 4 plots a sequence of running averages of the number of disks 
in the left chamber (each one taken over 1.25 X 10~^°sec), and provides an 
indication of the density fiuctuations for the chosen system parameters. 

Figure 5 shows the same quantities as figure 4, but examines them on 
a much longer time scale. The running averages of figure 5 are based on 
longer time intervals (25 X 10~^°sec), so the size of fiuctuations is accord- 
ingly reduced. The cumulative time average of the number of disks in the 
left chamber is a straight line that intersects the y-axis at the value of 0.68 
(normalized). This corresponds to 0.69 density (normalized) if we take into 
account that the unequal number of disks in each chamber changes the avail- 
able area. The simulation results are in good agreement with the theoretical 
estimate of 0.7 normalized density difference of section 3. 

5 The Loss of Ergodicity 

Our membrane appears to create a density difference without doing any work. 
In particular, it appears to establish the equivalent of a lower potential energy 
in the left chamber that creates a density difference in order to equalize the 
potential difference between the two sides. This can not be correct however. 
If we cut a hole in the membrane, a return flux of disks will result through 
the hole which can be used to convert thermal energy into useful work in 
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violation of the second law of thermodynamics. 

In order to rescue the second law of thermodynamics we can follow a 
number of different approaches. The first approach is to assume that the 
membrane can not be cut, that it must be a closed surface. In nature for 
example a chemical potential difference exists at the interface between two 
different materials such p-n semiconductors (see [3]). There is nothing in 
our model of the membrane however that implies that the membrane is an 
interface between two different materials and that it can not be cut. So 
we have to look elsewhere for an explanation of why the membrane can 
not violate the second law of thermodynamics. We start by examining the 
evolution of the disks in phase space. 

A system of N disks in two dimensions can be represented as a vector 
(. . . , Xi, Yi, Ui, Vi, . . .) of real numbers in R^^ phase space. This AN vector is 
the representative point of the system, and it specifies exactly the positions 
and velocities of the particles in the system at any given time (all particles 
have equal mass). As the system evolves in time, the representative point 
moves inside a constant energy subsurface of the R^^ phase space because 
the total energy is conserved. The total linear momentum is not conserved 
because wall collisions refiect a disk's momentum U- = —Ui or V/ = 
— Vi and because membrane interactions permute and/or refiect a disk's 
momentum according to equation 1. The total linear momentum is only 
conserved in a time average sense, and this has been checked by computer 
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simulations. 

The membrane and disks system is not ergodic because it does not spend 
equal times in equal regions of the accessible phase space (for example 
see [4, p.68]). If the representative point of the system visited regions of 
that have more disks in the left chamber as often as regions of that 
have more disks in the right chamber, then the time average density would 
be equal in the two chambers. Instead, the membrane and disks system 
approaches irreversibly a state of 0.7/0.3 density difference independent of 
initial conditions. 

In general there are many ways in which a system can fail to be ergodic. 
A trivial way to lose ergodicity in the context of billiard balls is to remove 
all interactions between the disks. Then the disks can not see each other 
and bounce between the walls of the container undisturbed, and the system 
does not attain a Maxwellian velocity distribution. This is a trivial example 
that shows that collisions between disks are necessary for modeling ideal gas. 
Binary collisions give rise to chaotic dynamics and allow a system of hard 
disks to explore fully its phase space. 

Binary collisions are not enough to guarantee ergodicity however. A sys- 
tem based on binary collisions and some other dynamics can fail to be ergodic 
if the additional dynamics introduce an attractor in phase space. In partic- 
ular this can occur when the evolution map Af : — )• does not conserve 
measure in phase space; that is the Jacobian determinant of the evolution 
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map is not equal to one (absolute value). The simplest example of a map 
that does not conserve measure is a two-to-one map, which means that two 
distinct representative points in are mapped onto the same point. In phys- 
ical space it means that two distinct trajectories of disks are mapped onto 
the same trajectory. An example is the one-way valve Maxwell's demon 
described in section 9. 

It turns out that the membrane and disks system also compresses phase 
space volume. The following three sections analyze this compression and 
demonstrate that the second of thermodynamics requires the incompressibil- 
ity of microscopic dynamics or an appropriate energy cost for compressible 
microscopic dynamics. 

6 Compression of Phase Space Volume 

To examine the compressibility of dynamics, we consider a membrane system 
that contains only one disk for simplicity. In other words, we have the same 
container and membrane as shown in figure 2 and we have a single disk inside. 
The phase space is three dimensional X, where X^Y is the position 
of the disk and is the velocity angle. We examine the compressibility of 
dynamics using this one disk system. 

It is easy to see that during free motion and collisions with the wall, the 
evolution of the system conserves phase space volume. The question is what 
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happens during the membrane-disk interaction. Referring to figure 3 we see 
that the transformation of the velocity angle of an incident disk is a linear 
map that has slope minus one, which suggests that phase space volume is 
conserved. This is incorrect however. In order to calculate the compressibility 
of dynamics we must consider the evolution of all three X, dimensions 
of the phase space together, and not only the dimension that is shown in 
figure 3. 

For concreteness we look at the time evolution of an infinitesimal phase 
space volume c<;(Xi, Ki, 0i) centered at the representative point Xi,yi,0i. 
We assume that the membrane is located at X = 0, and that Xi is near the 
membrane with Xi > , and that 0i is inside the interval 37r/4 < 0i < tt. 
We assume that after a time interval of 1.0 (in appropriate units) every 
point X, y, in the volume io{Xi^ Ki, 0i) has penetrated the membrane from 
the right and moved to the left of the membrane. We denote by X', Y\ 0' 
the image of X, under the evolution map, and we have the following 
equation, 

X' = X + cos0At, -sin0(l- At,) (7) 
where At, is the time it takes for point X, to move to the membrane 
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and is equal to —Xj cos 0. Therefore, 

X' = -sin0 -Xsin0/cos0 

Y' = y-X(l + sin0/cos0) -COS0 (8) 
0' = 37r/2-0 

To check whether the evolution map compresses the volume io{Xi , Ki, 0i) 
we evaluate the Jacobian determinant of the above equations. We find, 

-sin /cos - cos - X/cos^ 
-1 — sin0/cos0 1 sin0 — X/cos^0 
0-1 



sin0/cos0 = tan 



(9) 

For points X, inside the volume io{Xi, Ki, 0i), the angle is inside the 
interval 37r/4 < < tt. Thus, the Jacobian determinant is always less than 
one. In other words the evolution map compresses phase space volume when 
a disk penetrates the membrane from right to left. 

Figure 6 shows geometrically how the compression of phase space volume 
occurs. The figure is confined to two dimensions for graphical reasons. The 
two rectangles shown in solid lines correspond to phase space volumes that 
are mapped onto each other under the evolution map — they correspond to 
io{X, y, 0) with the angle chosen constant for all points. The vertical line 
X = of figure 6 corresponds to the membrane. It is easy to see that the 
edges of the two rectangles (the original rectangle and its image under the 
evolution map) are equal between the two rectangles, but the corresponding 



16 

angles are not equal, and hence the areas ol the two rectangles are not equal. 
Phase space volume is compressed when penetrating the membrane from 
right to left and expanded when penetrating from left to right. 

7 Other Membrane Maps 

The significance of compressible dynamics can be appreciated if we attempt 
to find a new membrane map that would result in a density difference while 
preserving phase space volume. It turns out that this is not possible, at 
least for piecewise differentiable maps. To see this, we seek a map /(©) 
mapping the impact velocity angle to a new velocity angle /(©) so that 
the condition of incompressibility (Jacobian determinant unity) is satisfied. 
Repeating the above steps, equations 8 and 9 using /(©), we find 

/(0) = sin-^(±sin0 + C) (10) 

where C is an arbitrary constant. If C is zero, /(©) corresponds to a trans- 
parent membrane (i.e. no membrane at all) or a completely impenetrable 
membrane (i.e. a wall). If C is non-zero, /(©) corresponds to a membrane 
that maps velocity angles in a non-linear fashion. In analogy with the mem- 
brane of equation 1, we apply the non-linear map /(©) of equation 10 to 
a selected region of velocity angles and block the remaining angles. In this 
way we hope that the probabilities of penetrating the membrane from the 
left and from the right will be unequal (see section 2). A simple calculation 
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however shows that the probabihties of penetrating the membrane from the 
left and from the right will be equal for all possible choices of the constant 
C in equation 10. 

For example let us suppose that (tt, tt — J) is an interval of velocity angles 
that penetrate the membrane from the right side, and (/(vr — J), /(vr)) is the 
image of (7r,7r — d) under the membrane map /(©) = vr — sin~^(sin0 + C) 
where C is a positive constant C <= 1 — sin J. Also let us assume that the 
remaining velocity angles are blocked and do not penetrate the membrane. 
Then the probability of an individual molecule to penetrate the membrane 
from the right (see section 3) is, 

P(R^L) = r \cose\de , (11) 

J TT — d 

and the probability to penetrate from the left is, 

P(L^H) = / cosOde . 12 

An elementary integration gives 

P(^R^L) = sin(7r — d) = sin d , 
P(L^R) = sin /(tt - d) - sin /(vr) = 

(13) 

sin[sin~"^ [sin(7r — d) -\- C]] — sin[sin~"^ [sin tt + C]] = 
sin d . 

Hence the probabilities of penetrating the membrane from the right and 
from the left are equal to each other. Therefore a membrane map (piecewise 
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differentiable) that obeys incompressible dynamics can not create a density 
difference. 

The above discussion confirms that the membrane of equation 1 achieves 
a density difference by compressing the phase space of the disks. The den- 
sity difference achieved by the membrane can be brought in accordance with 
the second law of thermodynamics by assuming that the membrane does a 
certain amount of work in order to compress the phase space of the disks. In 
other words, the second law of thermodynamics requires the incompressibil- 
ity of microscopic dynamics or an appropriate energy cost for compressible 
microscopic dynamics. 

8 Microscopic Rectifiers 

Another system that obeys compressible microscopic dynamics with an en- 
ergy cost is a pump of disks consisting of a microscopic trapdoor and external 
cooling. Such a pump is often used to explain thermalization [2, 8, 9]. In 
this section we examine the pump in terms of phase space compression. 

A microscopic trapdoor without external cooling is an example of a system 
that is designed to extract work from the thermal motion of disks by rectifying 
spontaneous variations in density in a system of disks. Such a trapdoor is 
also called a microscopic rectifier^ and it does not succeed as explained in 
references [2, 8, 9] because the rectifying mechanism becomes thermalized 



19 

and starts moving randomly in every possible way. In order to succeed the 
rectilying mechanism must be kept at a lower temperature than the system 
of disks. 

Interestingly the prevention of thermalization by a cooling process com- 
presses phase space volume. To see how compression of phase space volume 
occurs when a microscopic rectifier is cooled, we examine the trapdoor sys- 
tem of reference [2]. For simplicity we consider only one disk at first. We 
denote by x, j/, u, v the coordinates and velocity of the disk, and by X, U the 
position and velocity of the trapdoor. We assume that the cooling operation 
of reference [2] occurs at time Tc for < Tc < 1 (in appropriate units of 
time), and we consider the evolution of the system between times and 1. 
The new state is given by the following equations, 

x' = x-\-uTc-\-r]u (1 — Tc) 
y' = y + vTc + f]v [1 - Tc) 
u' = J] u 

(14) 

v' = rjv 

X' = X + UTc + eU (l-Tc) 
U' = eU 

where < e < 1 is the cooling parameter (a fixed number), and rj is chosen 
so that the total energy of the system is conserved. If m, M are the masses 
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of the disk and the trapdoor, the parameter rj is given by the formula, 



After some algebra we find that the Jacobian determinant of the evolution 
equation 14 is equal to e. In other words, phase space volume is compressed 
by a factor of e during a cooling operation. A similar calculation for a system 
containing a large number of disks, for example N disks, gives the following 
Jacobian determinant. 



which reduces to e (1 + 7(1 — e^)) if we assume that the total energy of the 
disks is much larger than the energy of the trapdoor, and that 7 is the ratio 
of the energy of the trapdoor to the average energy of the disks. Further if 
the trapdoor is colder than the disks, the energy ratio 7 is a small number, 
and the Jacobian determinant is approximately equal to e. Hence, cooling 
the trapdoor is accompanied by compression of phase space volume. 

A cooled rectifier is simply a pump and is not a threat to second law of 
thermodynamics [2]. The external cooling does work on the system and is 
responsible for compressing the phase space of the trapdoor and disks. When 
the cooled rectifier is enclosed in a larger system that is closed and isolated 
(for example if hot and cold reservoirs of finite heat capacity are used to 
perform the cooling), then the extended phase space evolves incompressibly 




(15) 



Jacobian 



= e 




(16) 
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and heat can only be converted to work while the system is approaching 
equilibrium. The membrane and disks system can also be viewed as a pump 
if we assume that the membrane does work in order to compress the phase 
space of the disks. 

9 Maxwell's demon 

An alternative way of bringing the membrane in accordance with the second 
law of thermodynamics is to view the membrane as a Maxwell's demon that 
interacts with the disks according to information about each incoming disk. 
For this purpose we review first the traditional one-way valve Maxwell's 
demon. 

Maxwell's demon is an imaginary being (or device) that operates a micro- 
scopic door between two chambers containing disks [5, 2, 6, 7]. The simplest 
version of the demon opens the door when a disk is coming from the right, 
and closes the door when a disk is coming from the left. The demon's op- 
erations lead to a density difference between two chambers, and the density 
difference can be used to extract work from the thermal motion of the disks. 
If the demon could operate in a complete cycle dissipating less energy than 
the work that can be extracted after the demon has finished its operations, 
then the second law of thermodynamics would be violated. This conundrum 
has inspired a large volume of literature aimed at exorcising Maxwell's de- 
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mon [6]. 

A popular way of exorcising Maxwell's demon (see [10]) is to assume 
that the demon operates its trapdoor according to information about each 
incoming disk, and to observe that the demon must erase the information 
that it has about the present disk in order to process the next disk. The 
erasure of information is irreversible because the demon maps two distinct 
disk trajectories onto the same trajectory every time it interacts with a disk 
by opening and closing its trapdoor. In particular after an interaction it is 
impossible to distinguish whether the disk came from the opposite chamber 
or whether the disk bounced off the demon's door. 

It is postulated that the irreversible erasure of information must be ac- 
companied by a minimum amount of entropy production which is ln[2) in 
the case of the one-way valve Maxwell's demon [10]. Thus the demon is 
brought into accordance with the second law because the reduction of en- 
tropy achieved by the demon is counterbalanced by an equal amount of en- 
tropy production that is necessary to implement the demon's irreversible 
operations. 

The membrane of equation 1 can also be viewed as a Maxwell's demon 
by imagining that the membrane-disk interactions are the result of a demon 
playing tennis with the disks. The demon moves a tiny racket up or down so 
as to intercept each incoming disk at the membrane line. In addition, the 
demon orients the racket in one of three possible orientations -|-45, —45, and 
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+90 degrees, so as to reflect the incoming disk according to equation 1 (also 
see flgure 2). 

The tennis demon difi'ers from the one-way valve Maxwell's demon in that 
the racket of the tennis demon must be positioned along a continuum of loca- 
tions depending on where the disk intersects the membrane line. In contrast, 
the door of the one-way valve Maxwell's demon must be positioned in one of 
two locations, closed or open. Another difference between the tennis demon 
and the one-way valve Maxwell's demon is that the tennis demon discards 
no information about the disks: The disk trajectories evolve bijectively so 
that the state of a disk (position, velocity) after a racket collision determines 
completely the state of the disk before the racket collision. In contrast, the 
one-way valve Maxwell's demon maps two disk trajectories onto the same 
trajectory every time it interacts with a disk. 

How can we then exorcise the tennis demon using information ideas? In 
section 6 we saw that the evolution of the disks compresses phase space 
volume. Now we may notice that the compression of phase space volume 
occurs because the demon chooses the location of the racket as a function of 
the incoming disk's position. In contrast, a racket that is positioned at a fixed 
location does not cause any compression of phase space volume. Therefore 
the problem must lie in choosing the racket locations. 

Suppose we consider a tennis demon that can only position the racket at a 
discrete number of flxed locations along the membrane line. In particular we 
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assume that the possible racket locations are uniformly spaced, and that the 
distance between successive locations is equal to Ax. If a disk trajectory does 
not fall exactly on one of the racket locations, the demon picks the nearest 
location possible. With the help of a diagram (see figure 7) we can see that 
no matter how small (but non-zero) the spacing Ax is, there are always disk 
trajectories that are mapped on top of each other. The discrete tennis demon 
does not compress the phase space of the disks in the continuous manner of 
section 6, but it does compress the phase space of the disks because of the 
many-to-one map of disk trajectories. 

The discrete tennis demon operates with finite information about each 
incoming disk, but it erases information irreversibly. The irreversible many- 
to-one map of disk trajectories disappears only in the limit of Ax going to zero 
(the spacing between racket locations). In this limit the discrete compression 
of phase space volume is replaced by the continuous compression of phase 
space volume of section 6. Moreover the reversibility of trajectories that is 
achieved in the limit comes at the expense of requiring the demon to operate 
with infinite information. Because Maxwell's demon can only operate with 
finite information (we can think of it as a microscopic computer), it follows 
that the tennis demon can not imitate the membrane of equation 1 reversibly. 
A tennis demon can only approximate the membrane of equation 1 using 
irreversible operations. 



25 



References 

[1] An interesting discussion of compressible dynamics has recently appeared 
in a paper that arrives at similar conclusions as ours, Kechen Zhang and 
Kezhao Zhang, "Mechanical models of Maxwell's demon with noninvariant 
phase volume", Phys. Rev. A 46 (8) 4598—4605 (October 1992). 

[2] P. A. Skordos and W.H. Zurek, "Maxwell's demon, rectifiers, and the second 
law: Computer simulation of Smoluchowski's trapdoor". Am. J. Phys. 60 
(10) 876—882 (October 1992). 

[3] R.P. Feynman, R.B. Leighton, and M. Sands, The Feynman Lectures on 
physics - Vol.3, (Addison- Wesley Publication, Reading, Massachusetts 1963, 
1965), chapter III-14, pp 14-9 — 14-10. 

[4] H.R. Pitt, Measure and Integration for Use, (Clarendon Press, New York 
1985), pp. 63—69. 

[5] W. Ehrenberg, "Maxwell's Demon", Sci. Am. 217 103—110, (November 
1967). 

[6] H.S. Leff and A.F. Rex, "Resource Letter MD-1: Maxwell's demon". Am. J. 
Phys. 58 (3) 201—209, (March 1990) ; Maxwell's Demon: Entropy, Infor- 
mation, edited by H.S. Leff and A.F. Rex, (Hilger/Princeton U.P., Princeton 
1990). 

[7] J.C. Maxwell, Theory of Heat, (Longmans, Green, London 1871), pp. 308 — 
309. 



26 



[8] R.P. Feynman, R.B. Leighton, and M. Sands, The Feynman Lectures on 
physics - Vol.1, (Addison- Wesley Publication, Reading, Massachusetts 1963, 
1965), chapter 1-46. 

[9] M. Smoluchowski, "Experiment ell nachweisbare der ublichen Thermody- 
namik widersprechende Molekularphanomene" , Phys. Z. 13, 1069 — 1080 
(1912). 

[10] C.H. Bennett, "Demons, Engines, and the Second Law", Sci. Am. 257, 108 — 
116, (November 1987). 



27 



Figure 1: A system of disks and a membrane can be viewed as a Maxwell's 
demon (see text). The membrane, displayed as a dashed line, interacts with 
the disks according to a time-reversible and energy conserving rule that cre- 
ates a density difference between the two chambers. 

Figure 2: The membrane interaction is shown graphically. The vertical solid 
line denotes the membrane, and the other solid lines denote a division of 
the plane into octants. The dashed lines correspond to disk trajectories that 
penetrate the membrane. 

Figure 3: The membrane interaction is shown as a map of the velocity angle. 
The impact velocity angle is mapped to the new velocity angle 0'. Both 
angles range from — tt to tt. The dashed line segments correspond to trajec- 
tories that penetrate the membrane, while the solid line segments correspond 
to trajectories that are mirror-reflected back. 

Figure 4: The membrane creates a density difference between two chambers 
initialized to have equal density. The smooth curve plots the cumulative time 
average of the number of disks in the left chamber, which increases from an 
initial value of 0.5 to approximately 0.7. The noisy curve plots a sequence of 
running averages of the number of disks in the left chamber, each one taken 
over f.25 x fO~^°sec. 
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Figure 5: The same quantities as in figure 4 are examined on a much longer 
time scale. The running averages are based on time intervals of 25 X f 0~^°sec. 



Figure 6: The membrane compresses the phase space of the disks. A re- 
gion of phase space that corresponds to disks with identical velocity angle is 
compressed when the points penetrate the membrane from right to left. 



Figure 7: A finite demon can position its racket at a discrete number of 
locations that are uniformly spaced along the membrane line. The demon 
picks the nearest location possible to position its racket, but there are always 
trajectories that are mapped on top of each other. The y-axis corresponds 
to the membrane and the solid lines are disk trajectories. 



